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Abstract. A very low primordial deuterium abundance 
of D/H ~ 1.5 10~ 5 has recently been proposed by Molaro 
et al. in the Lyman limit system with log Am — 18.1 cm~ 2 
at z a = 3.514 towards the quasar APM 08279+5255. The 
D/H value was estimated through the standard Voigt fit- 
ting procedure utilizing a simple one-component model of 
the absorbing region. The authors assumed, however, that 
'a more complex structure for the hydrogen cloud with 
somewhat ad hoc components would allow a higher D/H'. 
We have investigated this system using our new Monte 
Carlo inversion procedure which allows us to recover not 
only the physical parameters but also the velocity and 
density distributions along the line of sight. The absorp- 
tion lines of Hi, Cn, ClV, Si III, and Siiv were analyzed 
simultaneously. The result obtained shows a considerably 
lower neutral hydrogen column density log Am — 15.7 
cm~ 2 . Hence, the measurement of the deuterium abun- 
dance in this system is rather uncertain. We find that 
the asymmetric blue wing of the hydrogen Lya absorp- 
tion is readily explained by Hi alone. Thus, up to now, 
deuterium was detected in only four QSO spectra (Q 1937- 
1009, Q 1009+2956, Q 0130-4021, and Q 1718+4807) and 
all of them are in concordance with D/H ~ 4 10 -5 . 
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1. Introduction 

Accurate measurements of the hydrogen isotopic ratio 
D/H at high redshifts may allow us to test experimen- 
tally the basis of the standard model of big bang nucle- 
osynthesis (BBN) - its homogeneity. A homogeneous BBN 
implies a uniform distribution of the D abundance among 
the absorbing systems with low metallicity since 'no real- 
istic astrophysical process other than the Big Bang could 
produce significant D' (Schramm 1998, p. 6). 

First high quality spectral data of QSOs suggested, 
however, a dispersion in D/H = Adi/Ahi (the ratio of 
the Dl and Hi column densities) of about one order of 
magnitude (for a summary, see Buries et al. 1999). This 
finding provoked a lively discussion in the literature on 
inhomogeneous models of BBN (Dolgov & Pagel 1999 and 
references cited therein). But later, it was shown that a 
single D/H value of about 4 10~ 5 is sufficient to describe all 
observations available up to now (Levshakov et al. 1998a, 
1998b, 1999b). 

A new result by Molaro et al. (1999, hereafter MBCV), 
if confirmed, could challenge again the uniformity of the 
D/H space distribution because a very low deuterium 
abundance of D/H ~ 1.510 -5 was found in an extremely 
low metallicity system at z a — 3.514 towards APM 
08279+5255. MBCV consider the derived D abundance as 
a lower limit because their analysis was based on a simpli- 
fied one-component model of the absorbing cloud which 
failed to fit the red wing of the Lya line. They note that 
the observed complex structure of C iv and Si IV implies 
the presence of more than one component. They further 
state that 'additional components are required to repro- 
duce the extra absorption on the red wing of Lya' which 
would decrease the Hi column density for the major com- 
ponent leading to a higher deuterium abundance. 

Following the MBCV suggestions, we have analyzed 
the absorption profiles of Hi, Cn, Civ, Sim, and SilV 
by using our new Monte Carlo inversion (MCI) algorithm 
(Levshakov et al. 1999c, hereafter LAK). In the present 
Letter we show that the actual neutral hydrogen column 
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Fig. 1. Observed and synthetic line profiles of the absorption 
system at z a = 3.514. Note different intensity scales in panels 
h and f. For details, see text 
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density may be a factor of ~ 250 lower than the value of 
MBCV, a fact preventing any accurate measurements of 
the deuterium abundance in the z a = 3.514 system. 

The main difference between MBCV's analysis and our 
approach is that the former was performed in the frame- 
work of the microturbulcnt approximation using the stan- 
dard Voigt fitting procedure, whereas the latter is based on 
a more general mesoturbulent model which describes the 
process of line formation in the absorbing clouds more ade- 
quately (Levshakov & Kegel 1997). In both cases the same 
spectra of APM 08279+5255 obtained with the Keck-I 
telescope and the HIRES spectrograph by Ellison ct al. 
(1999) are utilized. 

2. Model assumptions and results 

Our model supposes a continuous absorbing gas slab of a 
thickness L (presumably the outer region of a foreground 
galaxy). The absorber is assumed to exhibit a mixture 
of bulk motions such as infall and outflows, tidal flows 
etc. Then the gas motion along a given line of sight may 
be described by a fluctuating (random) velocity field in 
which the velocities in neighboring volume elements are 
correlated with each other (while in the standard micro- 
turbulent approximation completely uncorrelated veloci- 
ties are assumed). The gas is tenuous and optically thin 
in the Lyman continuum. We are considering a compress- 
ible gas, i.e. the total number density of hydrogen «h is 
also a random function of the space coordinate, x. Fol- 
lowing Donahue & Shull (1991) and assuming that the 
ionizing radiation field is constant, the ionization of dif- 
ferent elements can be described by one parameter only 
- the ionization parameter U oc 1/nn- Furthermore, for 
gas in thermal equilibrium, Donahue & Shull give an ex- 
plicit relation between U and the kinetic temperature T. 
The background ionizing spectrum in our model is taken 
in the form given by Mathews & Ferland (1987). 

To estimate physical parameters and appropriate dis- 
tributions of the velocity v(x) and the normalized den- 
sity y(x) = nu(x)/no , n being the mean hydrogen den- 
sity, we used the MCI procedure developed on the ba- 
sis of the reverse Monte Carlo technique (Levshakov et 
al. 1999a). In our computations, the continuous random 
functions v(x) and y(x) arc represented by their sampled 
values at equally spaced intervals Ax, i.e. by the vectors 
{v\, ...,Vk} and {y\, ...,yk} with k large enough to repre- 
sent the narrowest components of the complex spectral 
lines. For the ionization parameter as a function of x, we 
have U{x) = U n /y(x) 7 with U n being the reduced mean 
ionization parameter defined below. We fix z a = 3.51374 
(the value adopted by MBCV) as a more or less arbitrary 
reference velocity at which Vj = 0. 

Our aim is to fit the model spectra simultaneously 
to the observed hydrogen, carbon and silicon profiles. In 
this case the model requires the definition of a simula- 
tion box for the six parameters (see LAK, for details) : 
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the carbon and silicon abundances, Zq and Zsi, respec- 
tively, the rms velocity and density dispersions, er v and 
<7 y , respectively, the reduced total hydrogen column den- 
sity Nh [= 7Vh/(1 + Cy) 1 / 2 ] , and the reduced mean ion- 
ization parameter U [= Uo/(l + 0- 2 ) 1 / 2 ]. For the model 
parameters the following boundaries were adopted : Zq 
ranges from 1CT 6 to 41CT 4 , Z si from 1CT 6 to 310" 5 , a v 
from 25 to 80 km s" 1 , er y from 0.5 to 2.2, A> H horn 5 10 17 
to 8 10 19 cm" 2 , and U a ranges from 5 10 -4 to 5 10~ 2 . 

Having specified the parameter space, we minimize the 
X 2 value. The objective function includes the following 
portions of the absorption profiles (labeled by grey lines 
in Fig. 1) which, after preliminary analysis, were chosen 
as most appropriate to the MCI fitting : for Hi Lya Aw 
ranges from —110 to —47 km s" 1 (the blue wing) and 
from 84 to 119 km s" 1 (the red wing, see Fig. la), for C n 
from -26 to 33 km s" 1 (Fig. lh), for CivA1548 from -51 
to 79 km s" 1 (Fig. lg), for Sim from -28 to 60 km s" 1 
(Fig. le), and for SiivA1393 Av ranges from -82 to 70 
km s -1 (Fig. Id). The necessity to choose these portions 
instead of the all available profiles was caused by a few dis- 
crepancies in the data. Thus two small spikes in C IVA1550 
at Aw = 15 and 23 km s" 1 are not seen in its stronger 
blue counterpart CivA1548, two minima in SiivA1393 at 
Av = —21 and 17 km s _1 are invisible in the profile of the 
red component SiivA1402. But in Fig. 1 (panels f and c) 
the observed CivA1550 and, respectively, SiivA1402 are 
shown together with the model spectra computed with 
the parameters derived from Lya, CnA1334, CivA1548, 
SiniA1206, and SiivA1393 fitting to illustrate the consis- 
tency. For the same reason the Ly/3 model spectrum is 
shown in Fig. lb at the expected position. The additional 
absorption on the blueward and the redward sides of Ly/3 
in the observed spectrum may be caused by Lya and/or 
Ly/3 forest lines at different redshifts. In Fig. 1, all model 
spectra are drawn by continuous thin lines, whereas filled 
circles represent observations (normalized fluxes). 

The MCI is a stochastic optimization procedure and 
one does not know in advance if the global minimum of the 
objective function is reached in a single run. Therefore we 
executed several runs for the given data set starting every 
calculation from a random point in the simulation box 
and from completely random configurations of the velocity 
and density fields. The results for five runs are listed in 
Table 1. The first solution (upper row in Table 1) was used 
to calculate the model spectra shown in Fig. 1, and the 
corresponding distributions of v(x), y(x), and T(x) which 
are presented in Fig. 2. 

The median estimations of the parameters give TVh = 
5.9 10 18 cm" 2 , N m = 5.3 10 15 cm" 2 , U = 1.6 10~ 2 , 
ct v = 51 kms-\ cr y = 1.1, [C/H] = - 1.8, and [Si/H] 
= - 0.7. The results were obtained with k = 100 and the 
correlation coefficients f v = f y = 0.95. 

The MCI allowed us to fit precisely not only the blue 
wing of the saturated Lya line but the red one as well. 



Table 1. Cloud parameters derived from the Lya and metal 
profiles by the MCI procedure (the total hydrogen column den- 
sity iV^f in units of 10 18 cm" 2 , the neutral hydrogen column 
density iV^ 1 in units of 10 15 cm" 2 , the mean ionization pa- 
rameter C/°2 in units of 1CP 2 , the rms velocity dispersion a v 
in km s" 1 , the reduced ^Xhi P er degree of freedom, v = 42, 
for the Lya blue and red wings) 
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* based on solar abundances from Grevesse (1984); 
[X/H] = log(JV x /JV H ) - log(iVx/iVH)0 



In addition we have reached a reasonable concordance 
between hydrogen and metal absorption lines. Although 
our model spectra of metal lines are consistent with the 
observed profiles within the 3-er uncertainty range, we 
found difficulties with the C n line fitting. The C n model 
spectrum lies systematically over the observed intensities 
(Fig. lh). We note that the measurement of the Cn line 
depends sensitively on the exact definition of the contin- 
uum because the line is very weak (it was not included by 
Ellison ct al. into the z a = 3.514 system). 

Nevertheless, let us assume for a moment that MBCV's 
identification of C II is correct, and let us further compare 
in some more detail the micro- and mesoturbulent results. 
Using the Doppler parameters 6c n = 17.8 km s _1 and 
&Hi = 21.0 km s" 1 measured by MBCV, it is easy to es- 
timate the turbulent velocity cr Vjmic = frturb/VS = 12.4 
km s _1 and the kinetic temperature T m j c = 8200 K. On 
the other side, the MCI kinetic temperature distribution 
(Fig. 2) docs not reveal values lower than 2 10 4 K. Be- 
sides, T m ; c = 8200 K seems to be too low for the pho- 
toionization heating. The mean kinetic temperature from 
the mesoturbulent solution T moso = 25500 K is, probably, 
more realistic. The four times higher mesoturbulent value 
of the rms velocity dispersion is also supported by obser- 
vations. Indeed, large broadening may be caused by the 
lensing effects since this QSO exhibits two components of 
similar intensity separated by ~ 0".4 (Irvin et al. 1998). 
The low T m i c and av, m ic seem to be an indication that the 
identification of Cn at z a = 3.514 is not very certain. 

3. Conclusions 

We have shown that the determination of the D abundance 
in the z a = 3.514 system is strongly model dependent. The 
Lya profile can as well be modeled with a considerably 
lower value of Ah i as compared to the MBCV model if one 
accounts for spacial correlations in the large scale velocity 
and density fields. This implies that the identification of 
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Fig. 3. Density-weighted radial velocity distribution functions, 
p(v), for the total and neutral hydrogen, Cn, Civ, Sim, and 
Si IV as restored by the MCI procedure. For comparison, in the 
upper panel the short dashed line shows p(v) for the homoge- 
neous microturbulent model adopted by MBCV 



and [Si/H] ~ -0.7. A metallicity as high as [C/H] ~ -1 
and a silicon overabundance [Si/C] = 0.5 — 1 dex has been 
measured in the metal systems at z a <~ 4 towards Q 0000- 
2619 (Savaglio et al. 1997). A similar silicon overabun- 
dance has also been observed in halo (population II) stars. 
The standard interpretation of these observations is that 
the metal-poor gas was enriched by Type II supernova 
nucleosynthesis products (e.g. Henry & Worthey 1999). 
Thus, our measurements are in agreement with these re- 
sults. 

Hence, it may be concluded that up to now deuterium 
lines have been identified convincingly in only four QSO 
spectra. Our previous D/H measurements in Q 1937-1009, 
Q 1009+2956, and Q 1718+4807 and the last one by Kirk- 
man et al. (1999) in Q0130-4021 are consistent with a 
single value of D/H ~ 4 10~ 5 . 
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the D I line in this system cannot be confirmed or ruled out 
without additional observations of the higher order Lyman 
series lines in order to constrain the total neutral hydrogen 
column density and the velocity field configuration. 

For this particular case, the presence of the extra- 
absorption exactly at the expected deuterium position 
Av = —82 km s _1 becomes evident from the upper panel 
of Fig. 3 where dotted and solid line histograms show, re- 
spectively, the density weighted radial velocity distribu- 
tions p(v) for the total and neutral hydrogen, whereas 
the dashed curve illustrates p(v) as adopted by MBCV 
(a Gaussian with cr v , m ic = 12.4 km s _1 ). The radial veloc- 
ity distribution of H I shows a blue-side asymmetry which 
may just mimic the deuterium absorption. The main max- 
imum of pe_i(v) is shifted to Av ~ 30 km s _1 with respect 
to MBCV frame of reference. The distribution phi(^) ex- 
hibits also an extended red wing ranging up to Av ~ 90 
km s" 1 . This asymmetry is also clearly pronounced in the 
Sim profile (Fig. le). 

The restored density and velocity fields reveal a very 
complex structure which is manifested in non-Gaussian 
distributions as shown in Fig. 3 for the total hydrogen 
density as well as for the individual ions. It appears that 
C iv is a good tracer for the distribution of the total hy- 
drogen density. 

Since the Voigt profile fitting analysis may produce 
very misleading results when applied to the case where 
the deviation from Gaussianity in the velocity distribu- 
tion becomes significant, the extremely low metallicity 
with [C/H] ~ -4.0 and [Si/H] ~ -3.5 reported by MBCV 
may be caused by their simplified model. For example, our 
analysis yields considerably higher values of [C/H] ~ — 1.8 
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